Cosmic Mach Number: A Sensitive Probe for the Growth of Structure 
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We investigate the potential power of the Cosmic Mach Number (CMN), which is the ratio between 
the mean velocity and the velocity dispersion of galaxies as a function of cosmic scales, to constrain 
cosmologies. We first measure the CMN from 4 catalogs of galaxy peculiar velocity surveys at low 
redshift (z G [0.002,0.03]), and use them to contrast cosmological models. Overall, current data is 
consistent with the WMAP7 ACDM model. We find that the CMN is highly sensitive to the growth 
of structure on scales k G [0.01,0.1] h/Mpc in Fourier space. Therefore, modified gravity models, 
and models with massive neutrinos, in which the structure growth generically deviates from that of 
the ACDM model in a scale-dependent way, can be well differentiated from the ACDM model by 
using future CMN data. 



Introduction — The Cosmic Mach Number (hereafter 
CMN) can provide a robust measure of the shape and 
growth rate of the peculiar velocity power spectrum of 
the galaxies in the universe. One considers a region of 
a given size r in the universe, and compares the bulk 
motion of the sphere with the random velocities within 
that region. The bulk motion provides a measurement of 
the forces on the region from irregularities external to it, 
so it measures the amplitude of perturbations on scales 
much larger than the region, whereas the random motions 
within the comoving region reflect gravitational pertur- 
bations on scales smaller than r. Thus their ratio, M{r), 
depends on the shape of the perturbation spectrum while 
independent of its amplitude. 

The concept was introduced by Ostriker and Suto in 
1990 [1] as a cosmological metric that would be relatively 
independent of the "bias" of the test particles being ob- 
served and also relatively independent of the then quite 
uncertain perturbation amplitude. They concluded that 
although the existing data were poor, they gave estimates 
of the CMN that appeared to be inconsistent with the 
then popular CDM model with r2,„ = 1 but seemed to 
prefer the open universe model instead. In a certain sense, 
the application of this test, correctly predicted the cur- 
rently best validated cosmological models with a value of 
Q,m in the range of 0.2 — 0.3. 

Subsequent to the original paper, Strauss et al (1993) 
[2] found again that the CDM models with fi^ = 1 re- 
mained inconsistent with the better data they used, but 
some non-standard models passed the test (see also [3]). 
Nagamine et al (2001) [4] looked at ACDM models and 
found better agreement, but the then ACDM model with 
= 0.37 again produced too high values of M over the 
range of 3 — 40Mpc/h, whereas a model with ilm ^ 0.2 
(actually closer to WMAP 7-yr constraint [5]) was more 
consistent with the observations. In addition, there are 
various other papers discussing the issues related to the 
CMN, such as non-linear clustering properties of dark 



matter on the CMN measurement [6] , and constraints on 
the CMN from Sunyaev-Zeldovich effect [7]. 

This history leads us to re-examine the issue in light 
of the much better knowledge now available from both 
the new peculiar velocity data and the range of models 
remaining plausible given the current cosmological con- 
straints. In this paper, wc will develop a new statistical 
tool to measure the CMN from the peculiar velocity sur- 
veys, and investigate the power of the CMN to distin- 
guish various cosmological models, especially in the as- 
pects of differentiating the ACDM model from variance 
with non-trivial growth function provided by Modified 
Gravity (hereafter MG) models, and from models with 
massive neutrinos. 

Statistics of Cosmic Mach Number — In linear pertur- 
bation theory, the power spectrum of the velocity diver- 
gence (6* = V ■ v) is related to the power spectrum of den- 
sity fluctuations via Pee{k,z) = f'^{k, z)P{k, z), where 
f{k,z) = ~d\ii5/dln{l + z), and 6 is the density per- 
turbation of matter. Since the data in our application 
are at very low- redshift, we assume that they have the 
same redshift z = throughout and drop the z depen- 
dence for brevity The mean square velocity dispersion 
(cr^(r)^ and mean square bulk flow {V'^{r)'^ in a window 
of size r can be calculated as [1, 4] 
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{V'{r)) dkPee{k)W{kr), 

i'y'ir)) dkPeoik) [1 - W{kr)] , (1) 



where M^(a;) = [3(sina;— xcosa;)/x^]^ is a top-hat window 
function. Note that W{x) 1(2; < 1) and W drops to 
quickly when a; > 1. Thus W effectively changes the 
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FIG. 1. (a): The 1-D posterior distribution of the CMN M{r) from 4 different catalogs, (b); The comparison between the CMN 
data (3cr CL.) and theoretic prediction: blue data points are the CMN data from posterior distributions shown in panel (a); red 
data points are the data used in Ostriker and Suto 1990 [1]. The black line is the CMN prediction from WMAP 7-yr best-fit 
values [5]; the green dashed line is calculated by using (1990) 'popular' CDM cosmological parameters Us = fim ~ 1 and h — 0.5 
[1]; the brown dashed and purple dot-dashed lines are the CMN from f{R) model with Bo = 10"'' and from ACDM model 
with neutrino mass ^ = 0.6eV. 



integral limits of the above formula as 
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where C = i% 



dkPgg{k) is a normalization constant. 
We can see that the velocity dispersion (CT^(r)) can be 
well approximated as a flipped signed large scale bulk 
flow {V^{r)) with an offset C. 

The CMN on different scales of the patches is defined 

as 



M{r) ^ V{V'{r))/{a^{r)) 
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(3) 

Thus it basically measures the shape of Pgg by contrast- 
ing j dkPgg{k) on large, and small scales. Using Eqs (2) 
and (3), one can reconstruct Pgg{k) up to an overall con- 
stant A, 



Pgg{k = l/r) = A 



M{r) 



dM{r) 



[1 + M2(r)]2 dr 



(4) 



Directly measuring Pgg from peculiar velocity surveys is 
challenging, but the measurement of CMN is much easier, 
as we show later. Fortunately, CMN has the same infor- 
mation as Pgg, except for the overall amplitude, which 
can be easily determined by probes such as the Cosmic 
Microwave Background (CMB). 

For a particular catalog with N objects, the CMN M 
can be written as M = |u|/cr*, where u denotes the bulk 
flow velocity, which is a streaming motion of galaxies to- 
wards some particular direction, and cr* stands for the 
small scale velocity dispersion. Unfortunately, neither u 
nor a* is a direct observable. For each galaxy peculiar 
velocity catalog, what we observe is the line of sight ve- 
locity Sn with measurement error (t„ for the nth galaxy. 



Then one can construct a joint likelihood function for u 
and M by contrasting Sn with the line-of-sight projection 
of the bulk flow fn,iUi- The uncertainty in {Sn — fn^iUi) 
is simply (c^ 4-0-2)5^ where ct* = |u|/A/ is given by the 
deflnition of the CMN. Therefore, the likelihood function 
takes the form of, 
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(5) 



One can then marginalize over the 3D bulk flow vector 
u to obtain the distribution of M for each survey. 

We use four different catalogs of galaxy peculiar ve- 
locity surveys, namely, SBF [8], ENEAR[9], Type la Su- 
pernovae (SN) [10], and SFI-I— |-f [11], to constrain the 
CMN by using Eq. (5). To calculate the characteristic 
depth f of each catalog, we use error-weighted depth as 
^ = J2n'^nr„/Y^n'^n, whcre Wn = I / {cr'i + ul). We 
marginalize over the 'bulk flow' velocity u in the 4-D pa- 
rameter space and obtain the 1-D posterior distribution 
of M as shown in panel (a) of Fig. 1. In the panel (a), one 
can see that the distribution of the CMN is very Gaus- 
sian, and the width depends primarily on the number of 
data entries in each catalog. In addition, since each cata- 
log probes the CMN on various depths, different catalogs 
form a complimentary set of tests of cosmic structures. 

In panel (b) of Fig. 1, we put together the old (1990) 
and current CMN data with various predictions com- 
puted from theoretic models: The [Blue] triangle data 
points are the current CMN values computed from likeli- 
hood (5) by using the four-catalogs, which provide more 
and deeper data than those used by Ostriker and Suto 
[1] in 1990 ([Red] square data). 
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FIG. 2. (a): Forecast la errors of the CMN from 6dF sur- 
vey, in which we assume the data is around our local region 
of the universe. Beside ACDM model, we plot Bo = 10"* 
([Brown] dashed line) and 10"'^ ([Red] dashed line) for the 
f{R) model, and "^m^ = 0.6eV ([Purple] dot-dashed line) 
and 0.3eV ([Green] dashed line), (b): APge/Pee for f{R) 
(dashed line) and massive neutrino (dot-dashed line) models, 
(c): Same comparison as (b) but for growth function A///. 



The ACDM model with WMAP 7-yr best-fit parame- 
ters (0„ = 0.271, h = 0.704, = 0.967, fib = 0.0455, 
[Black] solid line) is mildly consistent with the current 
CMN data at 3a CL. In comparison, we overplot the 
theoretical M{r) ([Green] dashed line) by using 1990's 
'popular' CDM parameters. One can immediately under- 
stand the reason why Rcfs. [1, 2] claimed that there was 
a strong conflict between the data ([Red] points) and 
popular CDM model ([Green] dashed line). There was 
in fact no inconsistency between data and a model with 
fl « 0.25 and ris — I. However, with the up-to-date data 
and strongly constrained ACDM cosmology, one can see 
that the CMN data is consistent with the ACDM pre- 
diction out to scale around 50Mpc/h. In addition, we 
plot the CMN for f{R) theory with Bq = 10"-* as brown 
dashed line, and the ACDM model with massive neu- 
trino 'Y^m^ = 0.6eV as purple dot-dashed line. These 
variations of growth factor exhibit some substantial dif- 
ference from the standard ACDM model. Note that the 
small scale CMN can be contaminated by the nonlinear- 
ity of the growth. For instance, using the fitting formulae 
of non-linear Pge proposed in [12], we find that CMN 
is suppressed about 21% (18%) on scale of 15Mpc/h 
(20Mpc/h). Thus nonlinearities must be considered when 
using the CMN data. Fortunately, there is a neat way to 
account for the nonlinearity when fitting models to the 
CMN data. From Eq (2), we can see that the nonlinear- 
ity only contributes to the constant C, and Eq (3) says 
that C can be treated as a nuisance parameter to be 
marginalized over in the parameter fit. 

To investigate the prospective accuracy of the CMN 
achievable in future surveys, we perform a forecast for 
the on-going 6dF pecuhar velocity survey [13]. The 
redshift distribution of galaxies for the 6dF survey is 
ng{z) = ^z^e-^^/^f)", where Zp ~ 0.0446, 7 ~ 1.6154 
and A ~ 622978 [13]. It peaks around z w 0.05 and ex- 
tends till z w 0.15. We assume that the 6dF survey can 
measure the peculiar velocity for 12000 brightest early- 



type galaxies, which is roughly 10% of the total sample, 
and they are located at z < 0.05, corresponding to the 
depth r < 150Mpc/h. We further assume that the mea- 
surement error for line-of-sight velocity is around 20%, 
which is a typical error for the fundamental plane dis- 
tance measurement. We divide these data into different 
redshift bins, and in each shell (r, r + dr), we calculate 
the Fisher Matrix value for the CMN Fmm from Eq.(5) 
which leads to the forecasted error of the CMN as 



C M{r 



M{r) 



1 



{u{r)/M{r)f 



(6) 



where N{r) is the number of data points in the shell 
(r,r -|- dr), and u[r) is the average bulk fiow magni- 
tude on depth r. We plot these forecast data in the 
panel (a) of Fig. 2. Comparing with panel (b) of Fig. 
1, We find the full range of the CMN data on scales 
[10,150] (Mpc/h) from 6dF can improve the constraint 
on the variation of the scale-dependent growth factor 
significantly. We summarize the experimental conditions 
for future experiment that can sharpen the CMN test: 
(1) there should be considerably more galaxy samples 
(> 10"*) on scales [10, 150]Mpc/h; (2) the smaller the 
measurement error (t„ is, and the larger sky area it cov- 
ers, the better it can reduce the overall error of M . One 
should also notice that, when the CMN data is used 
to constrain cosmology, the cosmic variance (CV) due 
to the limited volume of the survey should be taken 
into account. The cosmic variance crcv('') can be esti- 



mated as 



<ycv(r) 
M(r) 



T.4 (r) 



A{t) 



cb (r) 
B{r) 



where A 



and B stand for {V^{r)) and (o'^(r)) respectively, as de- 
fined in Eq (1), and a Air) and a sir) can be calculated 

■J „^\.,,P0e(k)\l - W(kr)]dk, where Ak and V 



are the width of the k bin and the volume of the survey 
respectively, and W is the same window function as in 
Eq (1). We found that for the current velocity surveys 
of SBF, ENEAR, SN and SFI-H+j., CV is 21, 19, 5 and 
10% of the measurement error. For the 6dF survey, CV 
is smaller than 5% of the statistical errors in all the bins. 

A sensitive test of growth of structure — Since the 
CMN measures the shape of the peculiar velocity power 
spectrum Pgg by design, it is sensitive to any distor- 
tion of Pgg. In the ACDM model, the growth is scale- 
independent. However, in the modified gravity models, 
and models with massive neutrinos, the growth is generi- 
cally scale-dependent, thus Pgg for these models is a dis- 
torted version comparing with the ACDM model, making 
the CMN an excellent tool to distinguish these models 
from ACDM. 

Let us consider the f{R) model and ACDM model 
with massive neutrinos as examples. In a viable f{R) 
model, where the Einstein-Hilbert action is extended to 
be a general function of the Ricci scalar R, the effective 
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FIG. 3. The 68 and 95% CL contour plots for Bo, neutrino 
mass and Qrn derived from a combined dataset of WMAP, 
UNI0N2 and 6dF surveys using Fisher forecast. See the text 
for details. 

value of Newton's constant GcS has both time and scale 
dependence, namely, GcS = fJ-{a, k)G, where /i(a, k) = 
(1 + |A2fc2a4)/(l + A^fc^a^), and G is the Newton's con- 
stant in general relativity (GR) [14, 15]. The only free pa- 
rameter is , which quantifies the Compton wavelength 
of the scalar field fn = df /dR and characterizes the scale- 
dependent growth rate. It is more convenient to re-define 
a dimcnsionless Bq which is Bg = 2HqX^ /c^, and Bq = 
for GR [15]. 

In panels (b) and (c) of Fig. 2, we show the relative dif- 
ference in Pge and / with respective to that in the ACDM 
model for a f{R) model with Bq — 10~^ (dashed lines). 
The growth rate for f{R) model is enhanced on scales of 
k > 0.01 h/Mpc, thus the integration of Pee (fc) cumulates 
the deviation in fc— space and exhibits the substantial dif- 
ference of CMN between ACDM and f{R) model. As a 
result shown in panel (a) of Fig. 2, a 0.01% of Bq in f{R) 
model can produce a 20% suppression in the CMN, which 
is potentially observable by 6dF. 

Similarly, the CMN is sensitive to the neutrino mass. 
When neutrinos became non-relativistic and the Universe 
was deeply in the matter dominated era, neutrino ther- 
mal velocities damped out the perturbation under the 
characteristic scale k^r — O.OlSfim {^miy/leV)h/Mpc, 
suppressing the power spectrum on small scales. On 
scales greater than fcnr , neutrinos affect the overall expan- 
sion of the Universe and therefore shift the peak of power 
spectrum to larger scales. In panels (b) and (c) of Fig. 2, 
we can see that neutrino with mass of 0.6eV can suppress 
the power spectrum on scales oi k > O.Olh/Mpc signifi- 
cantly, which exactly falls in the detection window of the 
CMN. Therefore, the cumulative 'integral' of Pee{k) in 
CMN can manifest the neutrino free-streaming effect by 
enhancing its value on all scales. Comparing with panel 
(b) of Fig. 1 and panel (a) of Fig. 2, we find that future 
CMN data is potentially able to distinguish the nonzero 
neutrino mass. 

Note that the possible degeneracy among cosmological 
parameters might dilute the constraints. Therefore we 
employ a Fisher matrix forecast [16] to study the power of 
CMN on cosmological parameter constraints. The Fisher 



matrix for CMN is, 

^."-E^Cov-.[M(..).M(.,)l^ (7) 

where p denotes the cosmological parameter, and Gov 
is the data covariance matrix for CMN. We include 
the cosmic variance of CMN in our analysis. We make 
a forecast for two models: A = {So, CP}, B = 
{STOiy,CP} where CP is a set of basic cosmological pa- 
rameters: CP= {fi,„, (78, -ffo, "s}- We choose the best- 
fit model for WMAP7 as a fiducial model [5], namely, 
CP= {0.266,0.801,71,0.93}, and choose Bq = 10"^ and 
EjTii/ = 0.05 eV as fiducial values for models A and B 
respectively. We also marginalize over the nuisance pa- 
rameter C, whose fiducial value was evaluated using the 
ACDM model. We also combine the Fisher matrices for 
the current data of CMB (WMAP [5]) and SNe (UNI0N2 
[17]) to improve the constraint (see details for CMB and 
SNe Fisher matrices in [14, 18]). The result is shown in 
Fig 3. As we can see, the constraints on Bq and neu- 
trino mass are largely improved when the CMN data from 
6dF is combined. Quantitatively, the 95% CL. upper lim- 
its on Bq and Sm^ are Bq < 0.4 [15] and Em^ < 1.3 
eV [5] using WMAP7 combined with UNI0N2 data, and 
the constraints can be tightened to < 5 x 10"^ and 
Stoi, < 0.65 cV when the CMN data from 6dF is in- 
cluded. 

Conclusion — We provide a statistical tool to measure 
the CMN. and further demonstrate that it is a sensitive 
probe of the structure growth. By design, the CMN is 
immune to the uncertainty in the overall amplitude of 
the density perturbation, and to linear galaxy bias. Also 
it is highly sensitive to any scale-dependent distortion of 
Pee since any small distortion can be amplified via the 
integral effect. Therefore it is an excellent tool to test 
alternative theories of gravity, and models with non-zero 
neutrino mass. 

We first perform a likelihood analysis of the CMN 
from the current peculiar velocity field data, and further 
confront the WMAP7 ACDM model and 1990s popular 
CDM model with the CMN data from current and old 
observations. We confirm that the ACDM model with 
WMAP 7-yr parameters is consistent with current CMN 
data at 3a CL. level. Based on our forecast for 6dF, we 
find that the CMN can improve the constraints on the 
modified gravity parameter Bq by 4 orders of magnitude, 
and it can also tighten the present constraints on the 
neutrino mass. The CMN information from future sur- 
veys, such as the Square Kilometre Array (SKA) [19], will 
be more powerful to constrain cosmologies, especially for 
modified gravity models, and models with massive neu- 
trinos. 
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